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A new study is reported of a lattice QCD calculation of the K + — » 7r + 7r° decay amplitude with the Wilson quark 
, action in the quenched approximation at /3 = 6.1. The amplitude is extracted from the K — > nir Green function, 
and a conversion to the continuum value is made employing a recent one-loop calculation of chiral perturbation 
theory. The result is consistent with the experimental value if extrapolated to the chiral limit. 
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1. Introduction 

It has long been a problem that the AI = 3/2 
K + — > 7T + 7r° amplitude calculated in quenched 
lattice QCD is about a factor two too large com- 
pared to the experimental value[0. In this ar- 
ticle we report a new study of this problem, in- 
corporating various theoretical and technical ad- 
vances in recent years for analysis. In particular 
we discuss in detail how one-loop corrections of 
chiral perturbation theory (CHPT) recently cal- 
culated^ affect physical predictions for the decay 
amplitude from lattice QCD simulations. 

Our simulation is carried out in quenched lat- 
tice QCD employing the standard plaquette ac- 
tion for gluons at (3 = 6.1 and the Wilson ac- 
tion for quarks. Two lattice sizes, 24 3 x 64 and 
32 3 x 64, are employed. We take up, down and 
strange quarks to be degenerate, and make mea- 
surements at four values of the common hop- 
ping parameter, n = 0.1520, 0.1530, 0.1540 and 
0.1543, which correspond to m^/m p = 0.797, 
0.734, 0.586 and 0.515. 



2. Extraction of decay amplitude 

The 4-quark operator most relevant for the 
AI = 3/2 K — > 7T7T decay is Q + — [(sd) l{uu) l + 
(su)i(Md)_L]/2. We extract the decay amplitude 
from the 4-point correlation function defined by 
M Q = (0\W+W o Q + (t)W K \0) where W ,+,k are 
wall sources for 7r°, ir + and K + . In our calcu- 
lations on a lattice of a temporal size T = 64, 
the walls are placed at the time slices t K + = 4, 
i 7r+ = 59 and t^o — 60. The mesons are all cre- 
ated at rest, and the 4-quark operator Q + is pro- 
jected to zero spatial momentum. 

We can use two types of factors to remove the 
normalization factors in Mq. If we define 

M w = (0\W Tr°(t)\0)(0\W + TT + (t)\0)(0\K(t)W K \0), 
M P = (0|Jf(t)Wx|0)(0|W + Wb7r + (t)7r (i)|0), (1) 



Mq/M) 
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we find Rw 
/(7r|7r|0) 3 • exp(£ - t + )A and R P = Mq/M p 
- (Tr+ir°\Q + \K+}/(ir\ir\Q) 3 for t K < t < t+,t , 
where A = — 2m 7! is the mass shift of the 
2-pion state due to finite lattice size effects 

In Fig. | we plot <7r|7r|0) 3 -R w and <vr|7r|0) 3 ■ R P 
at k = 0.1530 as a function of time t of the weak 
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Figure 1. R w ■ <7r|7r|0> 3 and i? P (7r|7r|0) 3 at k = 
0.153. Open and filled circles refer to data for 24 3 
and 32 3 lattices. 

operator. We clearly observe a non- vanishing 
slope for R\y, while Rp exhibits a plateau as ex- 
pected. The decay amplitude can be obtained 
by fitting Rw to a single exponential or Rp to 
a constant. We find the results to be mutually 
consistent within the statistical error. 

The one-loop renormalization factor relating 
the operator Q+ on the lattice and that in the 
continuum was obtained in Refs. We set 

q* = l/a or 7r j a as the matching point and em- 
ploy tadpole-improved perturbation theory with 
the mean-field improved MS coupling constant 
Quark fields are normalized with the KLM fac- 
tor §, y/1 - 3k/Ak c . 

3. Result 

In earlier calculations tree-level chiral pertur- 
bation theory (CHPT) formula was used to ob- 
tain the physical amplitude A phys from the lat- 
tice counterpart A lat . The formula takes the 
form AP h v s = (m 2 K - ml)/(2M%) • A lat where 
rriK — 497MeV and m w = 136MeV are physi- 
cal masses, and M„ is the degenerate mass of K 
and 7r mesons on the lattice. Clearly the physical 
amplitude should be independent of the lattice 
mass if the matching formula is exact. 

In Fig. H we compare results for the physical de- 
cay amplitude from a previous work|| carried out 
at P = 6.0 on a 24 3 spatial lattice with those of 
our simulation at (3 — 6.1. Our values obtained 
with the conventional quark normalization \/2/t 
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Figure 2. Comparison of our results normalized 
by the experimental value obtained with tree-level 
CHPT relation for q* = ir/a at (3 = 6.1 with those 
of previous work|Q] at j3 = 6.0 (crosses). Open 
symbols refer to results with the traditional nor- 
malization and no tadpole improvement, and 
filled ones are those with the KLM normalization 
and tadpole improvement. 

(crosses and open squares) are consistent with 
their results, but are larger than the experimental 
result roughly by a factor two. 

Let us also note with our results (circles and 
squares) that (i) the KLM normalization and tad- 
pole improvement of the operator have a signifi- 
cant effect on the amplitude, (ii) there is a clear 
dependence of the amplitude on the lattice meson 
mass M w , and (iii) a significant finite-size effect 
is observed between the two lattice sizes. These 
features show that the tree-level CHPT is inade- 
quate to extract the physical amplitude from lat- 
tice results. 

In Fig. ^ we show how predictions for the phys- 
ical amplitude change if we apply the one-loop 
formula of CHPT || to our results. Here A q and 
A cont are the cutoffs of CHPT for quenched and 
full theory. In converting to physical values we ig- 
nore effects of 0(p 4 ) contact terms of the CHPT 
Lagrangian since their values are not well known. 
An interesting point is that a size dependence 
seen with the tree-level analysis in Fig. || is ab- 
sent. Another important point is that the magni- 
tude of the amplitude decreases by 30 — 40% over 
the range of meson mass covered in our simula- 
tion. The amplitude depends significantly on the 
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Figure 3. Decay amplitude normalized by ex- 
perimental value for q* — n/a obtained with one- 
loop CHPT for A cont = 770MeV. Circles and 
squares refer to data for 24 3 and 32 3 spatial sizes. 
Open symbols are for A 9 = 770MeV and filled 
ones for 1 GeV. 

quenched lattice cutoff A q , in particular toward 
larger values of M n . 

Our results in Fig. || show that a sizable meson 
mass dependence still remains. This may be at- 
tributed to 0(p 4 ) contact terms which were not 
taken into account. If we assume that the contact 
terms of the full theory are as small as suggested 
by the phenomenological analysis ||, we can es- 
timate the physical amplitude by a linear extrap- 
olation of M% to the chiral limit. In Table |l| we 
list results for several choices of the cutoff and 
the operator matching point q*. As is expected 
from Fig. ||J variation of results on the choice of 
the cutoff parameters is small, being of the order 
of 10%. Allowing for this uncertainty, the values 
in Table [l] are consistent with the experimental 
value of 10.4 x 10~ 3 GeV 3 . 

4. Conclusions 

The present study has shown that the quenched 
result for the K + — > 7r + 7r° decay amplitude 
agrees with experiment much better than pre- 
viously thought, especially if the chiral limit is 
taken for the lattice meson mass. Away from 
the chiral limit, the amplitude depends signifi- 
cantly on the cutoff scales of CHPT, however. 
This problem is largely ascribed to a mismatch of 
chiral logarithms between the quenched and full 
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C + (7r7r|Q + |if)(xl0- 3 GeV 3 ) 

A cont A« 24 3 32 3 

(GeV) (GeV) q* = \ g* = f g* = \ q* = jj 

0.77 0.77 9.3(19) 10.2(21) 8.9(17) 9.7(19) 

0.77 1.0 9.4(13) 10.3(14) 8.8(11) 9.6(12) 

1.0 0.77 10.3(21) 11.3(23) 9.8(19) 10.7(21) 

1.0 1.0 10.4(14) 11.4(15) 9.7(12) 10.6(13) 



Table 1 

Results of linear extrapolation of the decay am- 
plitude to M% = 0. Statistical and extrapolation 
errors are combined. The experimental value is 
10.4 x 10~ 3 GeV 3 . 

theories In order to obtain K + — ► 7r + 7r° de- 
cay amplitude free from this uncertainty, we need 
simulations in full QCD. 
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